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ABSTRACT 

DNA methylation is an essential mechanism con- 
trolling gene expression during differentiation and 
development. We investigated the epigenetic regu- 
lation of the nuclear-encoded, mitochondrial DNA 
(mtDNA) polymerase y catalytic subunit (PolgA) by 
examining the methylation status of a CpG island 
within exon 2 of PolgA. Bisulphite sequencing 
identified low methylation levels (<10%) within 
exon 2 of mouse oocytes, blastocysts and embry- 
onic stem cells (ESCs), while somatic tissues con- 
tained significantly higher levels (>40%). In contrast, 
induced pluripotent stem (iPS) cells and somatic 
nuclear transfer ESCs were hypermethylated 
(>20%), indicating abnormal epigenetic reprogram- 
ming. Real time PCR analysis of 5-methylcytosine 
(5mC) and 5-hydroxymethylcytosine (5hmC) immu- 
noprecipitated DNA suggests active DNA methyla- 
tion and demethylation within exon 2 of PolgA. 
Moreover, neural differentiation of ESCs promoted 
de novo methylation and demethylation at the 
exon 2 locus. Regression analysis demonstrates 
that cell-specific PolgA expression levels were 
negatively correlated with DNA methylation within 
exon 2 and mtDNA copy number. Finally, using chro- 
matin immunoprecipitation (ChIP) against RNA 
polymerase II (RNApll) phosphorylated on serine 2, 
we show increased DNA methylation levels 
are associated with reduced RNApll transcriptional 
elongation. This is the first study linking nuclear 
DNA epigenetic regulation with mtDNA regulation 
during differentiation and cell specialization. 



INTRODUCTION 

Mammalian mitochondrial DNA (mtDNA) is a circular 
double-stranded genome of ~ 16.6 Kb encoding 37 genes 

(I) . These genes encode 13 polypeptides of the electron 
transfer chain, 22 tRNAs and 2 rRNAs, which directly 
or indirectly contribute to the production of ATP 
through the process of oxidative phosphorylation 
(OXPHOS). The remaining OXPHOS-associated genes 
and all of the transcription and replication machinery 
are encoded by nuclear DNA (2). Mutations, deletions 
and insertions to mtDNA lead to decreased OXPHOS 
capacity and mitochondrial function resulting in severely 
debilitating and often lethal pathophysiology (3). 

The mitochondrial specific DNA polymerase y (POLG) 
is the only known DNA polymerase that localizes within 
mitochondria (4,5). The POLG holoenzyme is composed 
of one catalytic subunit, POLGA, which is responsible 
for proofreading and DNA repair activity (6,7), and two 
POLGB accessory subunits, which maintain enzymatic 
stability and efficiency (8). POLG is essential for 
mtDNA replication and maintenance. Pathological muta- 
tions to PolgA are characterized by mtDNA deletion or 
depletion, which lead to respiratory chain deficiencies or 
premature aging phenotypes (9,10). Furthermore, ablating 

(II) or reducing (11,12) PolgA expression levels signifi- 
cantly decreases mtDNA copy number. 

The expression of PolgA and the numbers of mtDNA 
copy are strictly controlled during development. In preim- 
plantation embryos, the first mtDNA replication event 
is initiated at the blastocyst stage (13) and coincides 
with the up-regulated expression of POLG (13,14). This 
is restricted to the trophectoderm, which forms the 
placenta while those cells of the inner cell mass, which 
give rise to the embryo and embryonic stem cells 
(ESCs), do not replicate their mtDNA. As differentiation 
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and development proceed, significant changes in mtDNA 
copy number are observed (15-19), so that high-energy 
consuming tissues, such as muscle and nerve cells, 
acquire sufficient copies of mtDNA to meet their specific 
demands of ATP. The regulation of mtDNA copy number 
and PolgA expression during development are critical, as 
PolgA homozygous knockout mice die at E7.5-8.5 while 
heterozygous knockouts exhibit mtDNA depletion-type 
syndromes (11). 

Although some evidence indicates that PolgA expres- 
sion may be regulated by DNA methylation (20) and 
that mtDNA is itself methylated (21), the epigenetic mech- 
anisms that control the expression of PolgA and regulate 
mtDNA copy number remain to be determined. A CpG 
island within exon 2 of PolgA, which is hypermethylated 
during spermatogenesis, has been identified, however, its 
effect on PolgA expression and mtDNA copy number 
were not described (20). It is well-documented that 
cytosine methylation within promoter regions associated 
with transcription start sites mediates gene repression by 
recruiting methyl-binding domain repressor protein 
complexes or inhibiting the binding of transcription 
factors (22,23). However, recent studies have identified 
the importance of DNA methylation within (intragenic) 
and between (intergenic) gene bodies (22,24). The reason 
for DNA methylation at these sites remains unknown. 
Nevertheless, intragenic methylation has been linked to 
reduced transcriptional elongation (25-27) and gene 
silencing (28), while others have shown positive correl- 
ations between intragenic methylation and gene expres- 
sion (29,30). If such a mechanism regulated PolgA 
expression, this would present a unique regulatory mech- 
anism whereby an evolutionary distinct genome of bacter- 
ial origin is epigenetically regulated by a nuclear-encoded 
gene (31). 

We have determined the DNA methylation status of the 
intragenic CpG island within exon 2 of PolgA in gametes, 
pluripotent ESCs, differentiating stem cells and somatic 
tissues, and whether this intragenic DNA methylation 
correlated with PolgA gene expression. Using bisulphite 
sequencing, restriction enzyme digest real time PCR 
(qAMP) and methylated DNA immunoprecipitation 
(MeDIP) with real time PCR, we have identified tissue 
specific DNA methylation profiles and de novo intragenic 
DNA methylation during differentiation. We demonstrate 
that these DNA methylation profiles inversely correlate 
with the levels of PolgA expression and that there is a 
linear relationship between steady state levels of PolgA 
expression and mtDNA copy number. 

MATERIALS AND METHODS 

Cell culture 

Murine (Mus musculus) D3 ESCs (ESD3); somatic cell 
nuclear transfer ESCs (NT-ES) (32); CC9.3.1 ESCs 
(CC9"""); CC9.3.1 ESCs harbouring divergent Mus 
spretus (CC9 sprems ) or Mus terricolor (Dunni; CC9 dumi ) 
mtDNA (33) and induced pluripotent stem (iPS) cells 
derived from Quackenbush Swiss x Rosa26 {iPS^ s/R26 ) 
(34) and C57B1/6 x 129j/SvE Nanog-GFP + (iPS NGFP ; 



Stemgent, USA) (35) were cultured on mitomycin C 
treated mouse embryonic fibroblasts (MEFs) in ESC 
media composed of DMEM supplemented with 15% 
fetal bovine serum (FBS; Invitrogen, Australia), 0.5% 
(v/v) Penicillin-Streptomycin (Invitrogen, Australia), 
2mM L-glutamine (Gibco, Australia), 1% (v/v) 
non-essential amino acids (Gibco), 0.1 mM 
,8-mercaptoethanol (Gibco), with 1000 U/ml (lOng/ml) 
leukemia inhibitory factor (LIF; Millipore, Australia) at 
37°C in 5% C0 2 , as previously described (15). MEFs were 
cultured in ESC media in the absence of LIF and on 
gelatin coated culture dishes (0.1%). 

Tissue samples 

All mice (C57BLxCBA-Fl) were housed at the Monash 
University Animal House, according to Australian and 
Monash University ethical standards. Experimental 
procedures were approved by the Monash Medical 
Centre Animal Ethics Committee and were conducted in 
accordance with the Australian National Health and 
Medical Research Council (NHMRC) guidelines. 
Animals were euthanized by cervical dislocation and 
samples (oocytes, liver, spleen, kidney, muscle, heart and 
brain) were frozen immediately at — 80°C until required. 
MEFs (129/Sv) were donated by the Australian Stem Cell 
Centre. 

Embryo in vitro culture 

8-week-old female Fl mice (C57BLxCBA) were 
superovulated by intraperitoneal injection of 5 IU of 
human choronic gonadotropin (Chorulon; Invervet, 
Australia) and were mated with 12-week-old males. 
Embryos were retrieved 24 h post mating and in vivo 
fertilized embryos were collected from the oviduct follow- 
ing cervical dislocation of the female mouse. Culture to 
the blastocyst stage was performed, as previously 
described (36), in KSOM (Millipore, Australia) in 5% 
C0 2 , at 37°C. 

ESC injection into blastocysts 

Approximately 5-10 iPS yvG77p cells were injected into day 
3.5 (post-mating) blastocysts. Blastocysts were cultured on 
mitomycin C treated MEFs in ESC media at 37°C in 5% 
CO2. Outgrowths were selected after 5-7 days and 
manually passaged five times selecting only GFP + 
colonies. Further passaging was performed enzymatically 
using TrypLE (Gibco, Australia). 

Spontaneous and neural differentiation 

Spontaneous differentiation of ESCs was induced by 
culturing approximately 450 cells in 20 \A droplets of 
ESC media in the absence of LIF on the base of an 
inverted Petri dish for 48 h. Embryoid bodies were 
cultured in suspension for a further 5 days and then 
plated onto culture dishes in ESC media in the absence 
of LIF for 2 days. 

Neural induction was performed, as previously 
described (37,38). Briefly, ESCs were cultured at a 
density of 1.5-2 x 10 4 cm~ 2 on gelatin treated (0.1% w/v) 
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culture dishes in neural stem cell media (NSC) composed 
of DMEM/F12 (Invitrogen) supplemented with 1% N2 
(Millipore), 2% B27 (Millipore), 0.5% (v/v) Penicillin- 
Streptomycin (Invitrogen) and lOng/ul basic fibroblast 
growth factor (bFGF; Millipore, Australia). After 
3 days, the cultures were dissociated using TrypLE 
(Gibco) and plated on laminin (Sigma) coated (1 ug/cm 2 ) 
100 mm culture dishes at a density of 50000 cell cm -2 in 
NSC media. For neuronal differentiation, NSC media was 
changed every 48 h for a further 4 days, at which point 
cells were cultured in neuronal induction media (NIM) 
composed of Neurobasal Medium (Invitrogen) supple- 
mented with 2% B27 (Millipore), 0.5% (v/v) Penicillin- 
Streptomycin (Invitrogen) and lOng/ul brain-derived 
neurotropic factor (BDNF; Millipore). 

NCAM + cell isolation 

Spontaneously differentiated cells were triturated into a 
single cell suspension using TrypLE (Gibco) and passed 
through a 70 um filter. Single cell suspensions were 
incubated with ice cold isolation buffer (phosphate- 
buffered saline containing 0.5% bovine serum albumin 
and 2mM EDTA) supplemented with anti-PSA-NCAM 
conjugated to microbeads (1 in 11 dilution: Miltenyl, 
Australia) for lOmin at 4°C. PSA-NCAM labelled cells 
were isolated by magnetic separation Miltenyi MACS 
Cell Separation Columns (Miltenyi), according to the 
manufacturers instructions. 



RNA and DNA extraction 

Total RNA and DNA were extracted, according to the 
manufacturer's instructions, using the RNeasy Mini Kit 
(Qiagen, Australia) and the DNeasy Blood & Tissue Kit 
(Qiagen), respectively. The RNA samples were then 
treated with DNase I (Qiagen) for 15-30min at room tem- 
perature while the DNA samples were then treated with 
RNase solution (Qiagen) and proteinase K for lOmin at 
55°C. 

Reverse transcription was performed, as previously 
described (16). Briefly, 1 jxg/ja.1 of RNA, 5ul of lOx 
RT-Buffer, 1 mM of each dNTPs, 0.25 \ig Oligo (dT) 
primer, 1 U/ul RiboSafe, 1 U/ul of BioScript and 
autoclaved H 2 0 up to 20 ul were incubated at 42° C for 
lh and 72°C for lOmin. 



Generation of real time PCR standards 

Real-time PCR standards were produced using PCR. 
Breifly, 200 ng of cDNA or DNA were amplified in 50 ul 
reactions using 5x PCR buffer (Bioline), 1.5 mM MgCl 2 
(Bioline), 100 uM dNTPs (Bioline), 0.5 uM for each of the 
forward and reverse primers. Reactions were performed 
on a MJ Research PTC-200 machine and PCR products 
were resolved on 2% agarose gels at 100 V for 1 h. Primer 
sequences are listed in Supplementary Table SI. A series 
of 10-fold dilutions (1 ng/ul to 1 x 10~ 9 ng/ul) of the target 
specific PCR product was generated. 



Real time PCR: gene expression and mtDNA 
copy number 

Real time PCR contained 2ul of template (cDNA 
or DNA), 10 ul of 2x SensiMix (Bioline, Australia), 
0.33 uM of forward and reverse primers and ultrapure 
ddH 2 0 (20 ul total volume). Reactions were performed 
in a Rotorgene-3000 real time PCR machine 
(Corbett Research, Cambridge, UK) under primer 
specific conditions (see Supplementary Table SI). The 
number of mtDNA copies/cells and relative expression 
levels of PolgA were quantified against known standards, 
using real time PCR, as previously described (15). We 
assessed the reliability of fi-Actin, Gapdh and 18S as 
housekeeping genes. The expression of fi-Actin and 
Gapdh varied between cell types and tissues, although 
18S expression remained consistent across all samples 
analysed. 

Bisulphite sequencing 

Extracted DNA ( 1 rig) was incubated with 3 M NaOH for 
lOmin at 37°C and bisulphite treated, according to the 
manufacturers instructions (Human Genetic Signatures, 
USA). Mil oocytes (50 samples) and blastocysts 
(10 samples) were lysed and mixed with 1 rig of Herring 
sperm DNA (Sigma, Australia) prior to NaOH treatment. 
The samples were amplified using nested PCR primers 
(Supplementary Table SI; PolgA-Bimeth and PolgA- 
BiNest): 95°C for 4min followed by 35 cycles of 95°C 
for 1 min, 53° C for 1 min, 72° C for 1 min; and final elong- 
ation of 72°C for 4 min. PCR products were cloned 
using a pCR®2.1 vector (Invitrogen, Australia) and 
DH5a™ competent cells (Invitrogen). Individual clones 
were sequenced, as described in (39), using a 16-capillary 
3130xi Genetic Analyzer (Applied Biosystems, Australia). 
Data were analysed using a BIQ Analyzer (http://biq- 
analyzer.bioinf.mpi-inf.mpg.de/). 

Restriction enzyme digest and qAMP 

The methylation status of the CpG island in the 5' region 
of the PolgA gene (NT_039428: chromosome 7, region 
19 607 339-19 607 508), which has a variable DNA methy- 
lation profile during spermatogenesis, was analysed by the 
qAMP method (20). Genomic DNA was either digested 
with methylation-sensitive restriction endonucleases 
(Hpall, Notl), the methylation-dependent endonuclease 
(McrBC) or remained undigested (no enzyme sham diges- 
tion). The DNA samples were amplified by real-time PCR 
using primers defining a region including restriction sites 
in the second exon of the gene, which is included within 
the CpG Island (Supplementary Table SI; qAMP-PolgA). 
Shifts in Ct value (ACt) between the digested and un- 
digested samples were used to calculate the percentage 
of methylation at the assayed restriction sites. The ACt 
values were measured in triplicate. For the methylation 
sensitive enzymes {Hpall and Notl), percentage methyla- 
tion = 100 (2~ ACt ), and for McrBC, percentage methyla- 
tion = 100 (l-2" ACt ). 
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Immunoprecipitation of methylated DNA 

Immunoprecipitation of methylated DNA was performed, 
as previously described (40). Briefly, 10 fig of purified 
DNA was sonicated to produce random fragments 
between 200-1 100 bp. Sonicated DNA (3ug) was 
denatured at 95°C for lOmin and immunoprecipitated in 
500 ul of IP buffer (lOmM Na-Phosphate, pH 7.0, 
140mM NaCl, 0.05% Triton X-100) containing 2ug of 
either anti-5-methylcytosine (5mC; Active Motif, USA) 
or anti-5-hydroxymethylcytosine (5hmC, Active Motif) 
for 2h at 4°C. Samples were then incubated with 20 ul 
of Protein G Dynabeads (Invitrogen) for 2h at 4 °C and 
finally washed twice with 700 ul IP buffer for lOmin. 
DNA was purified using 7ul of Proteinase K (lOmg/ml; 
Qiagen) at 50° C for 3h and then precipitated using 1 jul 
Glycogen (20ug/ul; Invitrogen) and 2x volumes of 100% 
ethanol at — 80°C for a minimum of 2h. The enrichment 
of the nuclear-encoded PolgA and the mtDNA-encoded 
CytB, tRNA/COXl genes and the D-Loop was determined 
using real time PCR and standardized against /3-Actin 
(Supplementary Table SI). 

Chromatin immunoprecipitation 

Chromatin immunoprecipitation (ChIP) was performed as 
previously described (41) with minor modifications. Cells 
and tissue samples were cross-linked with 1% formalde- 
hyde for lOmin and quenched in 125 mM glycine (Sigma) 
for 5min. Heart samples were homogenized, passed 
through 200 um filter and centrifuged at 2000 r. p.m. -1 
for 5min at 4°C. Cross-linked samples were lysed at 4°C 
for lOmin and sonicated to fragment chromatin to an 
average length of 200-800 bp. Chromatin from 1 x 10 6 
cells was immunoprecipitated overnight at 4°C by 
adding Protein G Dynabeads (Invitrogen, Melbourne, 
Australia) pre-incubated with the antibody of interest. 
Antibodies used in ChIP assays were: polyclonal Rabbit 
anti-CTCF (Cell signalling, MA, USA); 4 ug monoclonal 
anti-RNA polymerase II (RNApII) clone 4H8; and 4ug 
monoclonal Mouse anti-Ser2 RNApII (Abeam, 
Cambridge, UK). Cross-links were reversed by incubating 
samples with 200 mM NaCl and 10 ul proteinase K(>600 
mAU/ml) at 65°C for 16 h. DNA was purified using the 
Qiagen PCR purification kit, according to the manufac- 
turer's instructions. Purified immunoprecipitated DNA 
and input DNA (1:100) were analysed using real-time 
PCR in triplicate, as described above. Primers are listed 
in Supplementary Table SI. Differences between cell types 
were calculated as percentage input (100 x 2 <Adjusted input 
(ct) - ip (ct»^ Yh ese va i U es were further normalized against 
negative controls. Enrichment values were converted to 
enrichment levels relative to exon 2. 

Statistical analysis 

Statistical analyses were performed using GraphPad Prism 
5.01 (GraphPad Software, Inc., San Deigo, CA) and 
Microsoft Excel. Data are expressed as mean ± standard 
error of the mean (SEM). All real time PCR data were 
normalized to 18S. Correlations were determined using 
Pearson coefficient and statistical differences between 



groups were determined using One-way or Two-way 
ANOVA in combination with Bonferroni post-hoc tests. 

RESULTS 

Tissue-specific DNA methylation within exon 2 of PolgA 

To determine whether DNA methylation to a region 
downstream of the PolgA transcription start site (i.e. 
intragenic methylation) regulates PolgA mRNA levels, 
we first analysed the levels of DNA methylation within 
exon 2 of oocytes, embryos and tissues. Bisulphite 
sequencing was performed on the region 1446-1616 of 
exon 2 (NCBI: NC_000073.5; chromosome 7, region 
19607 339-19607508) relative to the transcription start 
site of PolgA (20). Low levels of DNA methylation were 
detected in in vitro matured metaphase II (Mil) oocytes 
(8.44% ± 2.3; Figure 1A); in vivo fertilized blastocysts 
(0.71% (0.71; Figure IB); and ESCs (ESD3; 6.81% 
(1.85; Figure 1C). 

We then determined whether low levels of DNA methy- 
lation are present in differentiated cells by analyszing 
mouse embryonic fibroblasts (MEFs) and adult tissue 
samples (Figure 1D-K). All somatic samples displayed 
elevated levels of DNA methylation in exon 2 compared 
to Mil oocytes, blastocysts and ESD3 cells (P < 0.001; 
Figure 1L). A variable pattern of DNA methylation 
across the PolgA locus was observed for MEFs 
(45.23% ± 8.17; Figure ID) and liver (53.97% (2.18; 
Figure IE). Although a base-specific pattern of DNA 
methylation was observed in liver, there were no 
similarities in the DNA methylation patterns between 
cell types (Supplementary Figure SI). The highest levels 
of DNA methylation were detected in spleen 
(82.59% ± 6.93; Figure IF and L), which was significantly 
different to all somatic samples except muscle (Figure 1L). 
Our findings demonstrate a specific region within exon 2 
of PolgA is methylated in a tissue specific manner and that 
low or absent DNA methylation is associated with 
gametes and pluripotent cells while hypermethylation is 
prevalent in somatic tissues. 

Methylation of exon 2 correlates with reduced 
steady state mRNA levels of PolgA 

Since intragenic DNA methylation has been associated 
with regulating gene expression (27,28), we analysed the 
expression levels of PolgA. The relative expression of 
PolgA was significantly lower in somatic tissues 
compared with ESD3 cells (P< 0.001; Figure 2A). The 
range of mRNA levels observed in somatic tissues was 
between 69.9 (kidney) and 97.8% (liver) lower than that 
of ESD3 cells (Figure 2A). PolgA expression in spleen, 
brain, muscle and kidney were all statistically similar, 
while the levels in liver were significantly lower 
(P<0.05; Figure 2A). Regression analysis by Pearson- 
correlation coefficient demonstrated an inverse relation- 
ship between relative expression of PolgA and the 
corresponding DNA methylation profiles (R 2 = 0.5619; 
P = 0.032; Figure 2B cf Figures 1L and 2A). 

We then investigated whether the relationship between 
intragenic DNA methylation and PolgA expression 
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Figure 1. DNA methylation analysis of PolgA exon 2 in gametes, embryos, pluripotent and somatic cells. Bisulphite sequencing of the 14 CpG 
methylation sites in (A) Mil oocytes, (B) blastocysts, (C) ESD3 cells, (D) MEF (129/Sv), (E) Liver, (F) Spleen, (G) Heart, (H) Muscle. (J) Kidney 
and (K) Brain samples. Open and filled circles represent unmethylated and methylated CpGs, respectively. (L) The percentage CpG methylation per 
sequence determined by bisulphite sequencing (% mean ± SEM). Significant differences between cell types are: **T< 0.001. 



influenced mtDNA copy number. Analysis of mtDNA 
copy number demonstrated that the spleen contained 
251.9 ± 32.65 copies per cell, while liver, heart and 
muscle contained significantly more, namely 
2987.3 ± 564.3, 3956.02 ± 279.33 and 3199.17 ± 447.61 
copies per cell, respectively (P < 0.001; Figure 2C). 
ESD3 cells contained fewer mtDNA copies compared 
with all somatic cells and tissues except spleen 
(P < 0.001; Figure 2C). Despite tissue specific mtDNA 
copy number (Figure 2C) and intragenic DNA methyla- 
tion of PolgA (Figure 2B), linear regression analysis 



(Figure 2D cf Figures 1L and 2C) did not provide a sig- 
nificant correlation between mtDNA copy number and 
DNA methylation (R 2 = 0.0312; P = 0.6755). 

Intragenic DNA methylation in ESCs harbouring 
divergent mtDNA haplotypes 

Since there are significant interactions between nuclear 
DNA and mtDNA (3,11), we tested whether a more di- 
vergent population of mtDNA would affect the levels of 
DNA methylation to PolgA. We first analysed the levels of 
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Figure 2. DNA Methylation of exon 2 correlates with reduced steady state mRNA levels of PolgA. (A) Real time PCR quantification of PolgA 
expression in cultured ESD3 and MEF cells, and in liver, spleen, heart, muscle, kidney and brain samples, expressed relative to ESD3. (B) The 
relationship between DNA methylation levels from Figure 1L and the corresponding PolgA expression was determined using Pearson correlation 
coefficient (R 2 ). (C) MtDNA copies/cell in cultured ESD3 and MEF cells, and liver, spleen, heart, muscle, kidney and brain, as determined by real 
time PCR. (D) The relationship between the levels of DNA methylation from Figure 1L and the corresponding mtDNA copies/cell was determined 
using Pearson correlation coefficient (R 2 ). Values represent mean ± SEM and significant differences between cell types are: */"<0.05; **P<0.0\ and 
***/>< 0.001. 



DNA methylation in Mus musculus (CC9'""' V ) ESCs 
(Figure 3A and G) and found these levels were compar- 
able to ESD3 cells (Figure 3G cf Figure 1L). However, 
when this nuclear genotype harboured the more divergent 
Mus terricolor (Dunni) mtDNA (CC9 <A """) (33), which 
represents a divergence of 4.0 million years before 
present and a 1 59 amino acid difference (42), it displayed 
slightly lower levels of DNA methylation (4.16% ± 1.38; 
Figure 3B and G) compared with CC9"""' ESCs 
(9.89 ± 2.62; Figure 3A and G). 

The influence of nuclear reprogramming on 
intragenic DNA methylation 

Somatic cells can be reprogrammed to pluripotency 
by somatic cell nuclear transfer (32) and 
transcription-factor transduction, otherwise known as 
induced pluripotency (43). However, these approaches 
do not completely recapitulate the epigenetic and tran- 
scriptional signatures associated with ESCs (44) and, con- 
sequently, fail to correctly regulate mtDNA copy number 
during differentiation (16,31)- Therefore, we determined 
how somatic cell reprogramming affects DNA methyla- 
tion of PolgA in exon 2. Analysis of two independent 
induced pluripotent stem (iPS) cells lines, iPS QS (34) and 
■pgAOT/' demonstrated statistically elevated levels of 
DNA methylation for iPS QS cells (58.14 ±1.95; Figure 3C 
and G) compared with pluripotent ESD3 and CC9'"" S cells 
CP < 0.001). However, DNA methylation in iPS NGFP cells 
(20.13 ± 6.07; Figure 3D and G) was only elevated 



compared with ESD3 cells (P<0.05). Somatic cell 
nuclear transfer derived ESCs (NT-ES; 22.07 ± 5.14; 
Figure 3E and G) (32) were dissimilar to both ESD3 
and CC9 mus cells CP < 0.05; Figure 3G cf Figure 1L) sug- 
gesting that the exon 2 locus of PolgA is abnormally 
regulated in reprogrammed cells. 

We then asked whether the introduction of 
reprogrammed cells into the pluripotent environment of 
the inner cell mass of a blastocyst-stage embryo followed 
by their derivation into ESCs would modulate the DNA 
methylation signatures in exon 2 of PolgA. We micro- 
injected approximately 5 to 10 iPS^* 7 ^ cells into day 3.5 
blastocysts and the resultant Nanog-GFP + colonies were 
cultured for several passages. These iPS p ' Inj pluripo- 
tent stem cells displayed similar DNA methylation 
profiles (15.93 ± 4.51; Figure 3F) to non-injected cells, 
thus further demonstrating incomplete reprogramming 
(Figure 3G). Despite these differences, there were no 
CpG site-specific DNA methylation patterns identifiable 
between any of the divergent and reprogrammed cell types 
(Supplementary Figure S2A). 

The effects of DNA methylation on PolgA expression in 
divergent and reprogrammed cells 

CC9' m ' v , CC9 rf "" m ' and ESD3 cells (Figure 3H cf Figure 2A) 
all contained similar levels of PolgA mRNA levels while 
iPS Qs (p<o.05) and NT-ES cells CP < 0.01) had reduced 
levels. iPS /VGf7> cells also displayed lower levels of PolgA 
when compared with ESD3 cells (P < 0.05; Figure 3H). 
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We then collectively analysed the data presented 
in Figures 1L, 2A, 3G and 3H and confirmed the 
negative correlation (R 2 = 0.4597; P = 0.0055) between 
relative PolgA expression and DNA methylation (Figure 
3J). Taken together, DNA hypomethylation of the PolgA 
exon 2 loci is a common feature of pluripotent stem cells. 



However, the introduction of divergent Mus terricolor 
mtDNA into ESCs had no effect on the intragenic DNA 
methylation of PolgA and therefore PolgA expression 
while reprogrammed somatic cells displayed altered 
DNA methylation and PolgA expression compared with 
ESCs. 
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The effects of DNA methylation on mtDNA copy number 
in divergent and reprogrammed cells 

cc9 dunm ce j ls contained f ewer 

copies of mtDNA per cell 
(384.43 ± 18.39) than the parental CC9"""' cells 
(463.49 ± 18.92; i><0.05; Figure 3K). Furthermore, 
mtDNA copy numbers for iPS NGFP (441.61 ± 5.69; 
Figure 3K) and iPS ArG/rp - / '"' (435.07 ± 10.54; Figure 3K) 
cells were similar to ESD3 and CC9 mus cells while iPS QS 
(250.66 ± 5.85; P< 0.001) and NT-ES (234.20 ± 5.73; 
P< 0.001) cells contained fewer copies than ESD3 and 
CC9 mus cells (Figure 3K cf Figure 2C). Furthermore, 
using linear regression analysis, we observed a negative 
correlation between mtDNA copy number and PolgA 
expression (R 2 = 0.377; P = 0.0195; Supplementary 
Figure S2B). 

DNA demethylation within exon 2 

Recent studies have indicated that, during demethylation 
of DNA, 5mC is modified to 5hmC (45,46). Since 
bisulphite sequencing does not distinguish between 5mC 
and 5hmC, we quantified immunoprecipitated 5mC and 
5hmC MeDIP using real time PCR. MeDIP analysis of 
divergent ESCs, CC9 spretus (2.0 million years before 
present and 38 amino acid differences to the parental 
cells) and cC9 d "" m cells, demonstrated enrichment levels 
of 5mC (Figure 4A) and 5hmC (Figure 4B) within exon 2 
of PolgA, similar to CC9'"" 5 cells (Figure 4A and B). Liver, 
spleen and brain were all equally enriched for 5mC, while 
heart was further enriched (P < 0.001; Figure 4A). 
Enrichment for 5hmC was higher in heart and brain 
compared with liver and spleen (P < 0.001; Figure 4B). 
The enrichment levels of 5mC and 5hmC in tissues 
samples were elevated compared with ESCs (P< 0.001; 
Figure 4A and B). These enrichment levels are comparable 
to the DNA methylation levels identified by bisulphite 
sequencing, (cf Figure 1L; Figure 3G; R 2 = 0.4866; 
P = 0.0055; Figure 4C). We further confirmed the 
negative relationship between DNA methylation and 
PolgA expression (R = 0.6041; P = 0.0011; Figure 4D). 

5mC and 5hmC are present at very low levels in mtDNA 

As it has been suggested that mtDNA can undergo DNA 
methylation (21), we asked whether DNA methylation of 
the mitochondrial genome was directly linked to DNA 
methylation of PolgA in exon 2. Analysis of the levels of 
enrichment for the D-loop, CytB and tRNA/COXl 
regions of mtDNA from MeDIP samples revealed low 
levels of 5mC and 5hmC, similar to those identified for 
PolgA in ESCs and significantly lower than those observed 
in somatic tissues (P < 0.001; Figure 4E). 

Intragenic PolgA methylation and demethylation during 
spontaneous ESC differentiation 

During differentiation, mtDNA copy number is strictly 
regulated so that cells undergoing differentiation can 
increase their mtDNA copy number to meet their 
specific requirements for the generation of ATP through 
OXPHOS (15,16,31). We, therefore, determined whether 
intragenic DNA methylation of PolgA varied during 



differentiation. ESD3 cells were differentiated as 
embryoid bodies (EB) for 7 days and the DNA methyla- 
tion status of PolgA was determined using methylation- 
sensitive (Hpall and NotI) and methylation-dependent 
(McrBC) restriction enzymes followed by qPCR (20). 
We confirmed that this quantitative method for DNA 
methylation correlated with methylation profiles 
determined by bisulphite sequencing (Supplementary 
Figure S3A and B), and also validated the relationship 
between PolgA expression and DNA methylation 
(Supplementary Figure S3C). Within 24 h of ESD3 differ- 
entiation, DNA methylation at Hpall and McrBC 
sites increased to 100 and ~70%, respectively, with 
Hpall methylation decreasing to basal levels on Day 2 
(Supplementary Figure S3D). During differentiation, 
up to Day 7, methylation at NotI sites increased to 
levels <50%, while McrBC methylation decreased (Sup- 
plementary Figure S3D). Despite these changes in DNA 
methylation at the Hpall, NotI and McrBC sites, the 
overall mean methylation status was similar during early 
differentiation (Supplementary Figure S3E). These obser- 
vations demonstrate that the exon 2 locus of PolgA 
undergoes de novo methylation and demethylation 
during early ESC differentiation. 

mtDNA copy number and PolgA expression during early 
neural differentiation 

Since spontaneous differentiation of ESCs generates 
numerous cell types, we investigated whether the regula- 
tion of mtDNA copy number is developmentally regulated 
during neural and neuronal differentiation. In the first 
instance, ESCs were differentiated as a monolayer 
culture under defined culture conditions (37,38) for 5 
days to induce neural determination (Figure 5A). At this 
time point, we also isolated NCAM + cells in order to 
enrich for neural precursor cells, or allowed the cells to 
continue to differentiate for a further 7 days, under 
distinct culture conditions (37,38) (Figure 5A). In 
differentiating ESD3 cells, the mean mtDNA copies per 
cell was reduced after 5 days of neural-induced differenti- 
ation (both NSC and NCAM+ cells) when compared with 
their undifferentiated counterparts (Figure 5B cf 
Figure 2C; P< 0.001). However, as cells differentiated 
into neuronal cell types, mtDNA copy number increased 
to comparable levels detected in undifferentiated ESD3 
cells (Figure 5B cf Figure 2C). A similar pattern was 
also observed for the differentiation of CC9"'" S cells, 
although only ~35% reduction in mtDNA copies was 
observed in both NSC and NCAM + cells (Figure 5B cf 
Figure 3K; / > <0.01). In neural differentiation of ESCs 
harbouring Mus spretus mtDNA (CC9 spre,us ), mtDNA 
copy number also decreased compared with their undiffer- 
entiated counterparts and remained at lower levels at all 
time points of differentiation analysed (Figure 5B cf 
Figure 3K; P<0.01). Moreover, the above decrease in 
mtDNA copy number associated with early neural differ- 
entiation was not observed in mtDNA divergent 
NSC-CC9''"'"" cells (Figure 5B cf Figure 3K), while 
NCAM + -CC9 r/l """' cells isolated at the same time point 
contained significantly more mtDNA copies per cell 
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(Figure 5B; P < 0.05). These results indicate that the re- 
duction of mtDNA copy number during early neural dif- 
ferentiation is developmentally determined and increasing 
the genetic divergence between nuclear DNA and mtDNA 
alters this process. 

Similar levels of PolgA expression were detected in 
ESD3 and CC9' lwmi cells after 5 days of neural differenti- 
ation (NSC; Figure 5C). In NSC-CC9'" 1 " cells, the expres- 
sion levels of PolgA were significantly decreased compared 



with undifferentiated CC9'""' V ESCs (P< 0.001), while 
NSC-CC^"'""' cells expressed significantly increased 
levels compared with their undifferentiated counterparts 
(.P<0.01; Figure 5C cf Figure 3H). Nevertheless, the 
expression of PolgA was increased in NCAM + -CC9'"" V 
cells (P < 0.001) but reduced in both NCAM + -CC9 syr ''' MS 
(P< 0.001) and NCAM + -CC9' / ""'" (P< 0.001) cells 
compared with their non-isolated (NSC) counterparts 
(Figure 5C). 
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Neuronal differentiation over 12 days significantly 
reduced the expression of PolgA in ESD3 (P < 0.001), 
CC9' ms CP < 0.05), cC9 spretus (P<0.01) and cC9 dunni 
(P < 0.05) cells when compared with their undifferentiated 
counterparts (Figure 5C cf Figure 2A and 3H). DNA 



methylation levels at the exon 2 locus varied over the 
course of neural differentiation for all cell lines (Figure 
5D). Nevertheless, CC9 spre '"' s and CC9 dum,i cells displayed 
similar patterns of DNA methylation during neural differ- 
entiation while, for CC9 mus , a different pattern was 
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observed (Figure 5D). Since the transcript levels of PolgA 
are associated with the DN A methylation status of exon 2 
of PolgA in pluripotent and somatic tissues (Figures 2B 
and 3J), we again tested this relationship during 
differentiation. Regression analysis confirmed that DNA 
methylation within exon 2 of PolgA correlates with PolgA 
expression (R 2 = 0.4103; P = 0.0075; Supplementary 
Figure S4A; cf Figure 5C-D). Taken together, the 
changes in intragenic DNA methylation within PolgA 
during development provide further evidence of the epi- 
genetic regulation of PolgA expression. 

Intragenic DNA methylation of PolgA is correlated with 
expression levels in pluripotent, differentiating and 
somatic cells 

In order to further validate the association between PolgA 
expression and DNA methylation status of exon 2, we 
tested this relationship for all cells and tissues previously 
analysed. Regression analysis confirmed that DNA 
methylation within exon 2 of PolgA correlates with tran- 
script levels (R 2 = 0.4552; P = 0.0002; Figure 5E). 
Furthermore, using linear regression analysis, we 
observed a negative correlation between mtDNA copy 
number and PolgA expression (R 2 = 0.2244; P = 0.0145: 
Figure 5F; straight line). However, since the relationship 
between mtDNA copy number and PolgA expression also 
appeared exponential (Figure 5F; dashed line), we trans- 
formed the data and confirmed a significant relationship 
between mtDNA copy number and PolgA expression 
(R 2 = 0.5411; P = 0.0093; Supplementary Figure S4B). 

DNA methylation at the exon 2 loci is associated with 
reduced RNApII transcriptional elongation 

Having demonstrated the association of DNA methyla- 
tion with PolgA expression and mtDNA copy number, 
we further explored this relationship by examining tran- 
scription factor binding sites within exon 2 of PolgA. 
Using the JASPAR CORE transcription factor database 
(http://jaspar.genereg.net/cgi-bin/jaspar_db.pl) (47-50), 
we identified an evolutionary conserved putative CTCF 
binding site within this region (Supplementary Figure 
S5A). Since CTCF regulates transcription (41,51) and 
DNA methylation inhibits its binding (41), we determined 
whether CTCF binding regulates PolgA expression levels. 
Using the ChIP assay, we examined the enrichment of 
CTCF in CC9""' S ESCs and MEFs cells at the exon 2 
loci (Supplementary Figure S5B). Both cell types showed 
similar levels of CTCF enrichment within exon 2 and at 
sites upstream and downstream (Supplementary Figure 
S5C). Since CC9 mus and MEFs cells display different 
DNA methylation profiles and PolgA expression signa- 
tures, these observations suggest that CTCF is not 
involved in regulating PolgA expression. 

DNA methylation within coding sequences has been 
shown to alter transcript levels by influencing the recruit- 
ment or the transcriptional activation of RNApII 
(25,26,41). In order to determine how DNA methylation 
within exon 2 of PolgA modulates RNApII function, we 
quantified the enrichment of RNApII and the elongating 
RNApII phosphorylated at serine 2 of the carboxy- 



terminal domain (RNApIIS2). RNApII occupancy in 
CC9 m " ! , MEF, NSC-CC9'" 1 "' and heart samples was 
analysed upstream of the transcription start site within 
exon 2, and two sites downstream (Figure 6A). 
In CC9""" cells, RNApII (P < 0.001) and RNApIIS2 
(-P<0.01) enrichment was significantly increased down- 
stream of exon 2 (Figure 6B), while in MEF cells no dif- 
ferences were observed upstream or downstream of exon 2 
(Figure 6C). Furthermore, in NSC-differentiated CC9 mui 
cells, there was an initial increase in RNApIIS2 enrich- 
ment at a proximal site downstream of exon 2 
(P < 0.05), although this returned to similar levels 
observed at exon 2 further downstream (Figure 6D). In 
heart tissue, there was no difference in enrichment for 
RNApII between exon 2 and downstream sites, whilst 
RNApIIS2 enrichment was decreased downstream of 
exon 2 (Figure 6E; P < 0.05). These results indicate that 
DNA methylation downstream of the PolgA transcription 
start site reduces transcription by RNApII at the elong- 
ation stage. 



DISCUSSION 

The data presented here demonstrate a negative correl- 
ation between intragenic DNA methylation and the levels 
of PolgA expression, which provided an epigenetic mech- 
anism for regulating mtDNA replication. Through 
sequencing of bisulphite converted DNA, qAMP and 
MeDIP analyses, we have identified tissue-specific DNA 
methylation patterns within the gene body of PolgA, 
which increase as cells differentiate into specialized cell 
types. Changes in site-specific methylation levels were 
observed at this locus during early spontaneous differenti- 
ation of ESCs and 5hmC was enriched in post-mitotic 
tissues, which further suggest that changes in DNA methy- 
lation levels at the exon 2 loci in differentiated cell types are 
a result of active DNA methylation and demethylation. 

Experimental manipulation of peroxisome proliferator- 
activated receptor gamma coactivator (PGC) \-a and -fi 
has shown their levels of expression to be linked to 
mtDNA copy number [reviewed in (52)]. In this 
instance, hypermethylation of the PGC-la promoter 
correlated with reduced PGC-la mRNA and mitochon- 
drial density in type 2 diabetic patients (53). Despite these 
observations, PGC-la is not a critical factor in mtDNA 
replication since knockout mouse embryos develop to 
term and display only reduced mitochondrial biogenesis 
(54). The expression of the mtDNA replication factors 
Tfam (55) and PolgA (11) are essential for embryogenesis 
and regulating mtDNA copy number. Analysis of the 
Tfam promoter has suggested that DNA methylation 
within a nuclear respiratory factor- 1 (NRF-1) binding 
site may silence Tfam expression and reduce mtDNA bio- 
genesis (56). Others have suggested DNA methylation 
regulates PolgA expression during mouse postnatal germ 
cell development (20). Although, these studies suggest 
mtDNA replication may be epigenetically regulated, no 
association between DNA methylation and the expression 
of these replication factors or mtDNA copy number has 
been previously described. 
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DNA methylation of intragenic regions (within gene 
bodies) has been associated with negative (25,26,28) and 
positive gene regulation (29,30). The inverse relationship 
between PolgA exon 2 methylation and PolgA mRNA 
levels suggests intragenic methylation at this site 
suppresses mRNA expression levels rather than com- 
pletely repressing transcription. DNA methylation within 
gene bodies has been associated with alternative pro- 
moters and splicing of pre-mRNAs (26,41,57). The 
highly conserved transcription factor, CTCF, has a 
variety of roles in gene regulation (51) including mRNA 
splicing (41) and regulating transcription (51). Despite 
identifying an evolutionary conserved CTCF binding site 
within exon 2, ChlP analysis suggests PolgA regulation is 
independent of CTCF binding. Therefore, we further 
investigated the effects of DNA methylation within exon 
2 of PolgA by analysing the enrichment of serine 2 
phosphorylated RNApll, which is a marker of transcrip- 
tional elongation. Analysis following ChlP demonstrated 
that methylation within exon 2 is associated with lower 
levels of RNApll elongation downstream of this CpG 
island. Although the universal epigenetic function of 
intragenic methylation remains unclear and has been 
thought to differ between genes and genetic loci (22,58), 
our data demonstrate reduced transcriptional elongation 
regulates PolgA expression levels. 

The epigenetic regulation of chromatin structure and 
gene expression is essential for mammalian development 
and differentiation (59,60). Global DNA methylation 
analysis shows mouse ESCs are ~60% methylated (61) 
though stem cell genes, such as Oct4 (62), remain 
unmethylated in their promoters until differentiation 
induces gene repression through hypermethylation. 
Similarly, we observed low levels of DNA methylation 
within exon 2 of PolgA in oocytes, blastocysts and 
ESCs. These observations in gametes and pluripotent 



cells support a developmental relationship between 
cellular plasticity and mtDNA regulation. Indeed, 
siRNA knockdown of PolgA in mouse ESCs promotes 
differentiation and the loss of pluripotency (15). 

We have demonstrated that decreased PolgA expression 
during the early stages of spontaneous and neural-directed 
differentiation corresponded to increased levels of DNA 
methylation. During neural differentiation of mouse (15) 
and human (63) ESCs, mtDNA copy number reduces, 
which allows progenitor cells to establish the 'mtDNA 
set point', from which the mitochondrial complement 
may be expanded in a cell specific manner (64). The 
'mtDNA set point' was established in divergent ESCs 
although, during downstream differentiation, there were 
cell line specific differences in mtDNA copy number, 
PolgA expression and DNA methylation profiles for 
PolgA when compared with non-divergent controls. 
Thus, altering the mtDNA content of ESCs impacts on 
the expression of factors involved in mtDNA replication 
and maintenance during differentiation, and highlights a 
pathway, whereby changes in mtDNA content are 
communicated to nuclear DNA. 

By analysing mtDNA divergent ESCs and repro- 
grammed somatic cells, we have further demonstrated 
the importance of the strict epigenetic regulation of 
PolgA in pluripotent stem cells and during differenti- 
ation. Although pluripotent mtDNA divergent ESCs dis- 
played similar DNA methylation and PolgA expression 
levels to control ESCs, reprogrammed somatic cells had 
elevated DNA methylation profiles and variable PolgA 
expression. This suggests that incomplete reprogramming 
failed to reinstate the appropriate levels of DNA methy- 
lation to exon 2 of PolgA. We have previously 
demonstrated that reprogrammed somatic cells, including 
iPS and NT-ES cells, fail to regulate the levels of mtDNA 
copy number characteristic of early pluripotent stem cell 
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undergoing differentiation due to incomplete reprog- 
ramming (16). Indeed, abnormal epigenetic modification 
of chromatin structure is characteristic of somatic cell 
reprogramming (44,65) and is especially prevalent in 
NT embryos. This results in their failure to repress the 
expression of PolgA and they thus retain the ability to 
drive mtDNA replication during preimplantation devel- 
opment (39). 

As cell differentiation progresses, mtDNA copy number 
is regulated in a coordinated manner resulting in a termin- 
ally differentiated cell possessing the required number of 
mtDNA copies that reflects its cellular function (66). 
Through cumulative analysis of mtDNA copy number 
and PolgA expression levels, we ascertained a negative 
correlation linking the epigenetic regulation of PolgA to 
mtDNA copy number. MtDNA copy number was higher 
in somatic tissues than in ESCs but the expression of 
PolgA was higher in ESCs. The high DNA methylation 
and the enrichment of 5hmC in somatic tissues correlated 
with the reduced expression of PolgA indicating that 
demethylation may take place so that post-mitotic cells 
retain sufficient capacity to replicate mtDNA, as and 
when necessary. The higher levels observed in ESCs 
reflect their fast population doublings where dilution of 
mtDNA resulting from cytokinesis requires continual 
replacement. 

Recent data suggest that the mitochondrial genome is 
methylated (21) at levels higher than previously observed 
(67). DNA methylation is prevalent in prokaryotes (68) 
and eukaryotes (69) and is observed in chloroplast DNA 
(70). We observe very low levels of 5mC and 5hmC 
enrichment within mtDNA compared with the levels 
observed for PolgA in somatic tissues. DNA methylation 
of chloroplast DNA is not inhibitory to gene expression 
(70) while, in prokaryotes, DNA methylation is generally 
thought to limit the integration of foreign DNA (71), 
though, in some instances, it is associated with DNA 
mismatch repair, DNA replication and gene expression 
(72,73). It has been suggested that 5mC and 5hmC enrich- 
ment within mtDNA may influence mitochondrial 
transcription, although this association requires substan- 
tiation. Nevertheless, our data suggest that the epigenetic 
control of mtDNA replication during neural differenti- 
ation is coordinated by DNA methylation to exon 2 of 
the nuclear-encoded PolgA rather than by direct methyla- 
tion of mtDNA. 

In conclusion, we show that PolgA expression is nega- 
tively correlated with intragenic DNA methylation and 
reduced transcriptional elongation. Furthermore, we 
show that mtDNA copy number is regulated by PolgA 
expression in a cell-specific manner. This demonstrates 
that, in mammalian cells, the replication of the bacterial- 
originating mitochondrial genome is regulated by DNA 
methylation of a separate mammalian nuclear genome. 
These finding are especially significant in the context of 
nuclear-cytoplasmic coordination, whereby somatic 
reprogramming does not recapitulate ESC-like character- 
istics and divergent mtDNA can influence nuclear gene 
function during differentiation. 
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